Bare LO-Phonon Peak in THz-Emission Signals: a Dielectric-Function Analysis 
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We present a normal- mode analysis of coupled photocarrier-phonon dynamics in Te. We consider 
a dielectric function which accounts for LO phonons and the electron-hole gas within the Debye- 
Hiickel model and RPA. Our main finding is the existence of a bare LO phonon mode in the system 
even at high carrier density. This oscillation is an unscreened L-(q) mode arising from ineffective 
screening at large wave vectors. This mode is consistent with the bare LO-phonon peak in recent 
THz-emission spectra of Te. 

PACS numbers: 



The recent progress in ultrafast spectroscopy allows 
the investigation of photocajrier-phonon coupled THz 
dynamics in semiconductors^!. Ultrashort pulses enable 
the excitation and detection of coherent infrared lattice 
vibrations which is a source of THz radiation. THz emis- 
sion has been measured recently in semiconductors via a 
dipolercHitenna with detectable frequencies up to around 
5 THzl-a. Tellurium is a suitable semiconductor for these 
experiments because it possesses infrared-active modes 
within the antenna spectral range. 

The physics of THz emission in Te is straightforward. 
A laser pulse generates a non-equilibrium distribution of 
electron- hole pairs. The absorption length at the laser 
energies used in the experiments (~1.7 eV) is 40 nmEl. 
This small absorption length creates a highly inhomo- 
geneous longitudinal charge density and diffusion cur- 
rents arise. These currents lead to charge separation, 
due to the different electron and hole mobilities, which 
generates a longitudinal polarization field (Dember field). 
The Dember field drives LO phonons (Te is piezoelectric) 
whose macroscopic oscillatory polarization emits radia- 
tion. 

The emitted spectrum from the Te crystal consists of a 
broad peak at approximately 0.5 THz, followed by a dip 
to approximately zero at vto = 2-ArTHz and a second 
maximum close to vlo =2.82 THzLl,Q. The peak at 0.5 
THz_is related to the ultrafast buildup of the Dember 
fieldtj. The bare LO phonon mode is usually explained in 
terms of the lateral inhomogeneity of the excitation spot. 
In this picture the low-density "wings" of the spot gives 
rise to phonon-like long wavelength (q = 0) modes at the 
bare LO frequency [see L + branch in Fig. 1(a)]. 

Here we use a simple normal-mode analysis to gain 
further insight into the physics behind the bare LO peak 
observed in the emission spectrum of Te. Interestingly 
enough, we find that the longitudinal inhomogeneity of 
the electron-hole gas gives rise to bare LO phonon modes 
even at high carrier densities. This is in contrast with the 
usual picture involving lateral inhomogeneity of q = 
modes. We believe longitudinal inhomogenety (q 0) 
leading to ineffective screening of the phonon polariza- 
tion fields by the electron-hole gas provides an alternate 
(or complementary) description for the observed bare LO 



frequency. 

Dielectric function. In a first approximation the total 
dielectric function of the electron-hole-phonon system is 
given by 

e(q,u>) = £oo + /VX - V q P(q,u), (1) 
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with 



V q P{q,Lo) = ]T 
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where lo p i^ = J is the plasma frequency for elec- 

trons (i = e) or holes (i = h), m* is the effective mass 
( m */h = 0.067/0.53 mo), is the particle density, e is 
the electron charge, eoo = 10.9 is the dielectric constant 
for frequencies much above the reststrahl lo ^> loto = 
2ttvto, e o = 12-8 is the static dielectric constant, Tj is the 
carrier temperature and ks is the Boltzmann constant. 
Apart from a additive term, Eq. (|^) is essentially the 
Deby&JIuckel dielectric function within a hydrodynamic 
modelQ. Note that e(q, lj) takes into account the longitu- 
dinal inhomogeneity of the electron-hole gas through the 
q dependence in (|^) . Longitudinal inhomogeneity plays a 
fundamental role in the Te system because it gives rise to 
a driving Dember field which excites the infrared lattice 
modes. 

Normal modes. The zeros of the dielectric function 
(|l|) define the normal-mode frequencies of the system: 
e(q,u>) = uj(q,n), where n is the carrier density 
(we assume n = n e = n^). Figure 1(a) shows the usual 
dispersion relation for q = as a function of the carrier 
densitieso . The branches L + and L_ (solid lines) arise 
from the anti-crossing of LO phonons and plasmons. The 
important point here is that the L_ branch goes to the 
screened LO frequency vto as n increases. This is due to 
screening of the lattice polarization field by the plasma. 

Figure 1(b) shows the L_ branch for n — 10 18 cm~ 3 
as a function of q. For q — ► we recover the screened LO 
frequency since n is high. When q increases, however, 
screening becomes inefficient and the L_ branch shifts 
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to the bare LO frequency vlcM. Hence screening ineffec- 
tiveness is a possible mechanism to generate a bare LO 
peak in the spectrum of the electron-hole-phonon system. 

It is instructive to compare the Debye-Hiickel and RPA 
dielectric fwiction analyses. The RPA dielectric function 
is given byEI, 



V q P(q,to) 
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where is the sample volume, /i(£ik) is the Fermi distri- 
bution function for electrons or holes and £jk is the free 
electron or hole energy. Figure 2(b) clearly shows that 
RPA and Debye-Hiickel yield similar dispersion curves 
for L-{q). This is consistent with the high temperatures 
used in the calculation (T e/h = 7000/3600 K)M Note 
that we are not assuming q — > in the RPA contribution 

The main result of our simple normal-mode analysis 
- namely, that the longitudinal inhomogeneity of the 
electron-hole gas makes screening ineffective thus allow- 



ing for bare LO modes even at high densities "f-^lso follows 
from a more detailed numerical investigationEJ. We have 
solved a set of coupled differential equations consisting of 
Poisson's equation, a driven-harmonic-oscillator equation 
and transport equations which include inhomogeneity via 
diffusion terms. The numerical spectra are similar to the 
experimental ones, i.e., they also show a broad peak near 
0.5 THz, .-a dip at 2.6 THz and a second peak close to 
2.82 THzO. 

In conclusion, the dielectric function analysis presented 
here reveals that non-zero wave vector modes lead to bare 
LO phonon oscillations of the electron-hole-phonon sys- 
tem even at high carrier density. This is due to the sup- 
pression of the screening of the lattice polarization field 
by the inhomogeneous electron-hole gas. This screening 
ineffectiveness provides an alternative/complementary 
description for the observed bare LO frequency in THz 
emission experiments. 
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FIG. 1: Plasmon-phonon coupled modes for Tellurium as a 
function of carrier density [for q = 0] (a) and wave vector q 
(b). The lower branch L_ (a) tends to the TO frequency value 
for high densities due to the efficient screening of the lattice 
polarization field by the electron-hole gas. Ineffective screen- 
ing on the other hand, relevant for non-zero wave vectors q 
(b), makes the L_ branch move towards the LO frequency 
limit for increasing q's - even at high densities. This is in 
contrast with the q = mode in (a). 
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